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Double exchange mechanism is believed to favor transport along ferromagnetic directions, the 
failure of which in explaining the unusual resistivity anisotropy in iron pnictides is investigated. 
Several factors intrinsic to the microscopic mechanism of transport in iron pnictides are identified 
and analyzed, including the moderate Hund's coupling, low local moment, and presence of two 
anisotropic degenerate orbitals xz and yz. In particular, the substantial second neighbor hoppings 
are found to be decisive in giving results opposite to the double exchange picture. In high temper- 
ature nonmagnetic phase, orbital ordering is shown to give the right trend of resistivity anisotropy 
as observed experimentally, advocating its essential role in electron transport of iron pnictides. 

PACS numbers: 74.70.Xa. 74.25.F-, 75.30.Mb, 75.25.Dk 
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Close interplay among spin, orbital, and lattice de- 
grees of freedom together with the high-temperature su- 
perconductivity make iron pnictides a rich playground to 
study the effects of electron correlations. Unlike cuprates 
and manganites, all families of iron pnictides are metal- 
lic without doping. This greatly interests researchers 
in studying the transport property under different cir- 
cumstances and studying its relationship to all kinds of 
degrees of freedom [l|-l8|. Recently, an unusual resis- 
tivity anisotropy was observed in detwined underdoped 
Ba(Fei_a;Coaj2As2 [iHa] and was confirmed in other iron 
pnictides [71 l8(. The low temperature normal state of 
the material is an antiferromagnetic (AF) metal with 
coUinear (tt, 0) ordering. Experimentally, the resistivity 
was found to be larger along the ferromagnetic (F) di- 
rection than along the AF direction. This is thought to 
be in direct contradiction to at least two reasonings: (1) 
The ferromagnetic Hund's coupling between local mo- 
ments and the itinerant electrons indicates a double ex- 
change (DE) like correlation between spin configuration 
and transport, which implies a lower resistivity along the 
F direction [l|, |9|. (2) The scatte ring rate is expected to 
be larger along the AF direction [lO|, which is confirmed 
in optical conductivity experiment [J, |ll| . Another sur- 



prising result is that the resistivity anisotropy is largest 
close to the border between the AF phase and the super- 
conducting dome, where the magnetism and structural 
distortion are weaker [l|. 

Though conceptually puzzling, the early experimental 
results mentioned above are reproduced by several the- 
oretical calculations, including both mean field studies 
[id I12I [l3j and other more advanced techniques such as 
the dynamical mean field theory [l4l4l6| . Nevertheless, 
a clear physical picture concerning why the DE mecha- 
nism not working has not been addressed. Current un- 
derstanding of mean field results on resistivity anisotropy 
is only based on the topology and morphology of the 
Fermi surface [lO|, leaving the underlying mechanism 
untold. The failure of the interplay between ferromag- 
netism and conductivity in the AF metallic phase is thus 



hazy and still remains a puzzle. If the scattering rate 
is indeed stronger along the AF direction helping the 
electron transport along F direction [4!|, there must be 
a non-negligible competitor making the DE mechanism 
much less effective. Since several calculations via dif- 
ferent methods all observe this same unusual resistivity 
anisotropy, the underlying physics should be simple. It is 
thus interesting and timely to explore the applicability of 
the DE mechanism in iron pnictides, especially because 
the DE based model itself has been used successfully to 
unify the magnetic correlations in iron-based supercon- 
ductors [l7| and account for the novel spin dynamics [iSj . 
Besides, the reason of impairment of double exchange 
mechanism in transport may also happen in other mate- 
rials. In addition, orbital ordering is known to play an 
important role in many aspects of the system 19[422|. 
Surprisingly, mean field calculations have shown that or- 
bital ordering is irrelevant to 12| or even give opposite 
trend compared to [lO| the experimental results. It is ex- 
tremely desirable to reexamine the correlation between 
orbital ordering and the resistivity anisotropy presently. 

In this Letter, we report a theoretical study on the ap- 
plicability of DE mechanism in the electron transport of 
iron pnictides. A so-called orbital-degenerate DE model 
[iTl II8I is used to uncover the underlying mechanism of 
resistivity anisotropy both below and slightly above the 
structural and magnetic phase-transition temperatures. 
Although the DE effect is already weaken by the small 
Hund's coupling, low local moments and two degener- 
ate anisotropic orbitals, the relation between magnetism 
and transport anisotropy is still expected to follow the 
DE idea. The substantial second-neighbor hoppings are 
found to give rise to both a gapping feature in the band 
structure and a frustration in kinetic energy gain in trans- 
port along the F direction in the AF phase and is identi- 
fied as an important factor in competition with the DE 
mechanism. In the high-temperature nonmagnetic phase, 
the orbital ordering that could arise from a tiny exter- 
nal stress is shown to give the right trend of resistivity 
anisotropy and thus contributes positively to the experi- 



mental resistivity anisotropy. 

The model H=Ho+Hk+Hs+Hl contains four terms. 
Hq is the tight binding term up to second neighbor hop- 
Hfc is the local Hund's coupling between 
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localized spin and the itinerant electron spin, and is de- 
fined as Hk^-^KJ2^^^^^j^ Si • (4,M,a'^"/3'^i,M,/3)- A^ indi- 
cates the orbital index. K>0 is the strength of Hund's 
coupling, (Ti is the i-th Pauli's matrix acting in the spin 
subspace, and Si is the localized spin on site i. In order 
to maximize the DE mechanism, we take Si as classical 
spins [17[. Since we focus on the Q=(7r, 0) AF phase, we 
can take Si=S'zcos(Q • Ri). In the following, we rede- 
fine J=KS/2 to measure the strength of Hund's coupling. 
The first two terms of the model constitute the DE model 
in the limit of large Hund's coupling, which is well-known 
to facilitate transport along F directions |2J, [25 1. 

Hs is the superexchange term between localized spins, 
which competes with the first two terms to determine the 
ground state magnetic order and spin dynamics IT], ll8| . 
Since we concentrate on the electron transport in the 
fbced Q — (tt, 0) AF phase in this work, Hs is neglected 
in what follows. The final term is a crystal field splitting 
between the two orbitals as a result of the external stress 
field. It is written as H]^ — 5A^j(ni i — ni2), where 

ni,^—J2a 'H n o-'^i.A'.o- IS the number operator for electrons 
in orbital fj, and on site i. The lattice distortion also 
changes the hopping amplitudes along x and y directions 
[26| . However, we consider this as a minor effect in our 
present system of undoped or underdoped iron pnictides 
and neglect it in the following analyses. 

In the (tt, 0) AF phase, every unit cell contains two 
irons. Relabeling the electron annihilation operators on 
the two sublattices as ak,/i,cr and 6ic.p.<j. Take the a 



^k,2. 



..&L.J- 



spin basis as V'L={"k,i,<T' ^k,i, 

ing Hs , the original model is written as H=Hq 

HL=J2wa'^laK'>Pka, where 
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K(k) = 
-2i2COs(fcj^C2) - f- - (J J, /:^t(k) = -2^2 cos(fcj,C2) - 

Y + '^J: /3(k) = -2t2COs(fca;Ci)-4i3COs(fc:j;Ci)cOs(fcyC2), 
/4(k) = -2iiCOs(fc2;Ci)-4i3COs(fc:rCl)cOs(fcyC2), /5(k) = 

— 4i4sin(fc2,ci) sin(fcyC2). Ci (02) is lattice constant along 
X (y). Along {k^, 0) or (0, ky), /5(k)=0, the two eigen- 
energies for the two orbitals are written compactly as 

A^,± = -2t^cos(V2)-(-iry±V'-^^ + /^\2(k), (2) 

with /i = 1 (/i = 2) for the xz (yz) orbital. From the 
above expression, the two orbitals open gaps both along 
the kx and the ky directions. Along fc^;, the positions of 



gap opening of the two orbitals are far away from each 
other by a distance of |2(ti — ^2) — A|. For experimentally 
realistic J, no full gap opens along k^, which is consis- 
tent with a topological argument by Ran et al ;27]. This 
Dirac cone structure appears along k^ as an accidental 
crossing of the energy bands corresponding separately to 
the two orbitals. However, it could be proved that for 
A=0 and small J the gap positions for the two orbitals 
always coincide with each other along ky and are cen- 
tered around ^1^2/^3- This orbital-cooperated gap would 
highly enhance the corresponding resistivity. 

For the hopping integrals, we consider mainly the pa- 
rameters proposed by Raghu et al [23] and fix the largest 
hopping integral as ti—QA eV J7,, 19] . The ratios among 
the four hopping integrals ^1:^2:^3:^4 are 1.3:(-l):(-0.85):(- 
0.85). Results for another set of parameters used by Yin 
et al [i3|, in which ii:t2:i3:t4=4:1.3:(-2.5):0.7, are found 
to give qualitatively the same results. These two sets of 
parameters agree qualitatively in the sense that the near- 
est neighbor intra-orbital hopping integral of the xz {yz) 
orbital along the y (x) direction ti is larger than along 
the X {y) direction ^2, a peculiar feature resulting from 
the strong hybridization of the iron orbitals with the p 
orbitals of surrounding As ions [19| . 

The DE mechanism, first introduced for manganites, 
applies when the ferromagnetic Hund's coupling J is 
much larger than the hopping amplitudes and for a partly 
filled band [2J, |25|. Though the Hund's coupling for iron 
pnictides is not as large, it is not obvious why a result op- 
posite to expectation based on DE is observed. To have 
a feeling on the global behavior of the model, we show 
in Fig. 1(a) the resistivity ratio Py/ Px in a wide range 
of J, at 50 K in the low temperature AF phase. Two 
electrons per iron is considered for the undoped system 
|23l |28J, |29| . The finite temperature resistivity (conduc- 
tivity) is obtained from standard Kubo's formula [30|. 
For low J the resistivity anisotropy agrees with exper- 
iments. But when J increases to an unphysical value 
for iron pnictides, the resistivity anisotropy changes to 
the opposite behavior as compared to experiment, which 
is however in agreement with expectation from the DE 
mechanism. One peculiar feature on Fig. 1(a) is for very 
large J the resistivity anisotropy of the undoped system 
changes back to the behavior same as in experiment, a 
result specific to the undoped and extremely low doped 
systems and easily understood from the Pauli's exclusion 
principle for electrons in a half filled band. 

Before answering the question why undoped and 
slightly electron doped iron pnictides behave differently 
from DE predictions, we first calculate the resistivity 
anisotropy for a single orbital model with isotropic in- 
plane hoppings up to nearest neighbor, that is to set 
ti = i2 and ^3 = ^4 = 0. The results show that even 
for extremely small J, the resistivity is smaller a long the 
F direction, following exactly the DE prediction 2J, |25| . 
Thus, the small magnitude of J is not the crucial factor 
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FIG. 1: (a) Ratio of resistivity as a function of J between 
py along the F direction and px along the AF direction, for 
three typical dopings at 50K. (b) py/ Px for isotropic one or- 
bital (labeled as 'lorb') and two orbital (labeled as '2orb') 
models with 2NN hopping, at 50 K. (c) and (d) are disper- 
sions along kx and ky for J = 0.1 eV both in the presence 
(solid lines) and absence (dashed lines) of 2NN hoppings, re- 
spectively. The horizontal dotted lines in (c) and (d) mark 
the chemical potential. 



to bring about the observed anisotropy. 

However, the inclusion of substantial second nearest 
neighbor (2NN) hoppings in combination with the small 
magnitude of J literally induces the unusual resistivity 
anisotropy in iron pnictides, which is further enhanced by 
the large anisotropy in the intraorbital nearest neighbor 
(NN) hoppings. To see the effect of 2NN hoppings, we 
show in Fig. 1(b) the resistivity anisotropy for a model 
with one isotropic orbital (^2=^1=0. 4 eV, t3=-0.25 eV, 
t4=0) and two isotropic orbitals (^2=^1=0. 4 eV, i3=-0.25 
eV, ^4=0. 07 eV). Though with very different line shapes, 
the same qualitative behavior as in Fig. 1(a) is obtained. 
From Fig. 1(b), it is also obvious that a nonzero 2NN 
inter-orbital hopping t^ tends to reduce the magnitude 
of the resistivity anisotropy. 

To understand the role played by 2NN hoppings, first 
notice that in a square lattice the effective hopping along 
each of the two in-plane directions have contributions 
from both the NN and the 2NN hoppings. In the (tt, 0) 
AF phase, the NN and 2NN hoppings along x both con- 
nect AF bonds while they connect F and AF bonds re- 
spectively along y. When J is very large, the 2NN hop- 
pings together with the NN hoppings along the AF di- 
rection are effectively blocked, the DE mechanism simply 
prevails. However, in the region of small J, a more in- 
tricate competition between the maximization of kinetic 
energy gain and minimization of Hund's energy appears. 
Since hopping along the F bonds does not cost or gain 



magnetic energy, maximization of kinetic energy gain for 
J <C ii prefers a nearly non-polarized on-site spin dis- 
tribution. Whereas along the AF direction, the up and 
down spin electrons each experiences effectively a peri- 
odic shallow potential wells of depth 2 J, both of two-site 
period and shifted from each other by one site. This pair 
of periodic potential wells tends to drive the formation 
of notable local spin polarizations, which can lower the 
Hund's energy at the cost of kinetic energy loss. Since 
for small J the kinetic energy gain is dominant, for negli- 
gible 2NN hoppings the spin polarization follows largely 
that of the F direction. However, the additional mag- 
netic energy gain from hoppings along AF 2NN bonds 
enhances the spin polarization, as verified through ex- 
plicit calculations. This in turn frustrates hopping along 
the F bonds and finally impairs the DE for large enough 
is. We believe that the above is the physical picture why 
substantial 2NN hoppings invalidate the DE mechanism 
for very small J. 

The effect of the 2NN hoppings could be further un- 
derstood from the band structure. For clarity, the band 
structures for J = 0.1 eV in the presence and absence 
of the 2NN hoppings are shown together by solid and 
dashed lines in Figs. 1(c) (1(d)) along k^ {ky). In the 
absence of 2NN hoppings, the energy bands are more 
metallic along ky than along kx since no energy gap opens 
in the whole energy range in the former. When the 2NN 
hoppings are turned on, the band structures for undoped 
and slightly doped systems would be gapless along kx 
but would open a gap along ky as we have mentioned 
following Eq.(2). This peculiar change in the direction 
of gap opening tendency follows exactly the change from 
DE behavior to anti-DE behavior and thus provides a 
very clear picture how 2NN hoppings invalidate DE. The 
gap along ky exists once lia/iil > 1/2, providing a good 
reference critical value for the 2NN hoppings. 

Another important feature of iron pnictides different 
from standard DE systems (e.g., manganites) is that the 
two relevant orbitals are anisotropic with respect to x 
and y, which is manifested through |ti/t2| > 1- Com- 
paring Fig. 1(a) and Fig. 1(b), this intrinsic anisotropy 
enhances the resistivity anisotropy and gives rise to the 
peak structure in the small J region on Fig. 1(a). The 
low J peak is far less obvious in the absence of this 
anisotropy, as could be seen from Fig. 1(b). 

Former mean field calculations found smaller resis- 
tivity anisotropy for high moment states than for low 
moment states [lOJ, ll4|- According to Fig. 1(a), the 
reduction of resistivity anisotropy in the high moment 
phase indicates that iron pnictides should belong to the 
region changing from anti-DE behavior to DE behav- 
ior. We notice that there are iron based superconductors 
which are considered as having larger effective J, such as 
FcTci_a;Sea; [17|. One prediction worthy of future exper- 
imental inspection is that opposite resistivity anisotropy 
with larger resistivity along the AF direction in the un- 
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FIG. 2: (a) Ratio of resistivity between the direction with 
short lattice constant (py) and the direction with long lattice 
constant (px), for three typical dopings. Temperature is taken 
as 150K, without of magnetization, (b) The orbital polariza- 
tion S12 (in percentage) close to the chemical potential for a 
series of increasing energy windows SE. 



derdoped region of these materials might be observed. 

One common feature in iron pnictides is that a struc- 
tural transition always occurs at a higher or the same 
temperature as compared with the AF transition [1], in 
which the orbital ordering and fluctuations are consid- 
ered by many to be essential 19l42l|. It is thus very sur- 
prising that former mean-field studies showed that orbital 
ordering is anti-correlated with the resistivity anisotropy 
in the AF phase 10|, [12 1 . Experimentally, the same resis- 
tivity anisotropy is observed in the high temperature non- 
magnetic phase close to the AF and structural transition 
when an external stress field is applied to induce a lattice 
distortion following the structural anisotropy in the AF 
phase [l[ . Here we would like to single out the role played 
by the orbital ordering in this phase. In the low doping 
region, we expect crystal field effect to be dominant, and 
consider a dominance of the xz orbital for a lattice with 
a larger lattice constant along x, implying A < in our 
model [26J. As shown in Fig. 2(a), orbital ordering gives 
the right trend of resistivity anisotropy compared to ex- 
periment. We then calculate the orbital polarization for 
states close to the chemical potential, which are most 
relevant to electron transport, for a typical doping and 
a series of increasing energy windows SE in Fig. 2(b). 
Denoting the integrated density of states of xz (yz) or- 
bital within SE as Sni {Sn2), the orbital polarization S12 
is defined as {5ni — 5n2)/{5ni -f 5n2). It is clear that 
(5i2 close to the chemical potential has the same char- 
acter as the total orbital polarization set by the sign of 



A. Note that fluctuations via nematic correlations [1|.|15| 
neglected here might be important to provide a complete 
picture for the high temperature resistivity anisotropy 
and deserves future investigations. 

The above correlation between orbital polarization 
and resistivity anisotropy is encoded in the intrinsic 
anisotropy of the two degenerate orbitals. As could be 
seen from Eq.(l) for J = 0, the effective intraorbital hop- 
ping along kx is —2{t2 + 2t^) cos{kxc) for xz orbital and is 



parameter sets considered, \ti +2t2,\ < 1^2 + 2^3! 
So in a qualitative sense, a dominant xz polarization close 
to the chemical potential implies an easier hopping and 
conduction along x direction. In comparison, yz has a 
larger effective hopping along the y direction and thus 
facilitates the conduction along y direction. 

In conclusion, the puzzling interplay between magnetic 
order and electron transport in iron pnictides is iden- 
tified quite different from DE-dominated systems (e.g., 
manganites). The small Hund's coupling, low local mo- 
ments and intrinsic anisotropy in the NN hoppings all 
make the DE a bad description. In particular, the sub- 
stantial 2NN hoppings are found to be decisive in com- 
peting with the DE mechanism, greatly frustrating the 
kinetic energy gain in transport along the ferromagnetic 
direction. We have also found in the high temperature 
nonmagnetic phase that, an xz dominated polarization 
which might arise as a result of an external stress could 
give the same trend of resistivity anisotropy as observed 
in experiments, advocating the possible essential role of 
orbital ordering in the resistivity anisotropy. 

We reemphasize that 1^3/^1! ~ 0.5 is a good reference 
critical value of the 2NN intra-orbital hopping to im- 
pair the DE mechanism (for small J), the exact value 
of which depends on t^ and J. First-principles 2NN 
hoppings for all three t2g (two e^) orbitals are larger 
(smaller) than half of the corresponding NN hoppings 
[19l] . Since the Fermi surface of iron pnictides are dom- 
inated by t2g orbitals, a model containing more orbitals 
should not change our previous arguments. Similar im- 
pairment of the DE mechanism might be observed in 



other i2g-dominated systems [31|. SrCrOs is also an AF 



metal with (tt, tt, 0) magnetic order similar to iron pnic- 
tides ^32.] . However, our first-principles calculations give 
2NN hoppings smaller than half of NN hoppings for all 
five orbitals, predicting that DE should persist in that 
system. This work was supported by the NSC Grant No. 
98-2112-M-001-017-MY3. Part of the calculations was 
performed in the National Center for High-Performance 
Computing in Taiwan. 
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